On the phases of Polonium 
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The thermodynamical properties of the main phases of metallic polonium are examined using 
Density Functional Theory. The exceptional nature of the solid-solid phase transition of a to P 
Po is underlined: it induces a lowering in symmetry, from cubic to rhombohedral, with increasing 
temperature. This is explained as the result of a delicate balance between relativistic and entropic 
effects. Overall agreement with existing experimental data is very good by state-of-the-art standards. 
The phonons of Po present Kohn anomalies, and it is shown that the effect of spin-orbit interactions 
is the inverse of that in normal metals: due to the non-spherical nature of the Fermi Surface, 
spin-orbit effects reduce nesting and harden most phonon frequencies. 

PACS numbers: 61.66.Bi, 61.50.Ah, 63.20.dk, 64.70.kd 

A frantic race is underway to find exotic new phases 
of elemental solids under pressure [H 0, 0, 0). With the 
advent of extreme diamond-anvil cell experiments with 
pressures beyond 1 MBar, a whole slew of new phases 
can be synthesized: with increasing pressure, the distance 
between atoms decreases, electron hybridization changes, 
and the nature of chemical bonds as well, changing co- 
ordination numbers and bonding directionality, orbital 
occupations, and so on. This is particularly dramatic in 
the recent demonstration that simple metals can become 
semiconducting if valence and core orbitals are forced to 
hybridize, through the application of pressure 0, 

Comparatively much rarer are purely temperature in- 
duced solid to solid phase transitions. Although about 
half 5] of the elements present some form of solid-solid 
transition with T, most are in very limited pressure 
ranges, the critical pressure being more or less constant 
with T. Almost all of these transitions either increase or 
maintain the level of crystal symmetry (usually hexago- 
nal to cubic, or cubic to cubic). This is intuitive, as tem- 
perature will increase inter-atomic distances, and tend 
to rotationally average orbital configurations, leading to 
high-symmetry (most often BCC) then liquid structures. 
The exceptions are the usual suspects: elements with 
complex bonding (boron, sulfur, phosphorous and bis- 
muth), and a few of those with delicate electronic shell 
effects, such as Ce or Dy. In these systems, an increase 
in temperature can give rise to a solid phase with lower 
symmetry. 

Polonium also falls into this category. It is one of the 
strangest elements of the periodic table, and the only 
element which adopts the simple cubic (SC) [(| structure 
at ambient pressure (a-Po, joined by a few others such as 
Ca-III and As-II|7j at higher pressure). Upon increasing 
temperature, between 290 and 330 K,8j, Po transforms 
to the (3 phase, with a rhombohedral structure. 

Despite its intriguing nature, experiments on Po are 
very difficult and rare, and little is known about its prop- 
erties. Only a few laboratories in the world can synthe- 



size Po in sufficient quantities to perform characteriza- 
tion. Due to their high radioactivity (with a half life of 
at most 100 years, depending on the isotope), samples 
heat and decay very quickly into Pb or Bi. 

By chemical analogy, Po should behave like Se and Te, 
its 4p and 5p cousins, and adopt a trigonal spiral struc- 
ture at low T. However, it was shown recentlvjij, [Io| that 
the SC structure is stabilized over the spiral structure by 
relativistic effects (though there is debate over whether 
the scalar relativistic (SR) or spin-orbit (SO) terms dom- 
inate - see comments to Ref. [9J) . It was predicted that a 
simple rhombohedral structure would be stabilized under 
pressure. A qualitative explanation for the appearance 
of the SC structure is that the valence p states do not 
hybridize with s, because of strong relativistic repulsion 
of the s 11 , 12j . The p states can therefore sustain purely 
right-angle bonding and an SC structure. Something sim- 
ilar probably happens in Ca and As, due to the relative 
(but not necessarily relativistic) changes in the s and p 
energy levels under pressure. 

In this Letter, the vibrational and thermodynamical 
properties of crystalline Po are calculated from first prin- 
ciples. First, the effect of the spin-orbit interaction on the 
phonons of a-Po is shown to be counter-intuitive, harden- 
ing phonons and reducing the nesting of the Fermi surface 
(FS), the opposite behavior of simple metals. Thermody- 
namical quantities are calculated for the low temperature 
phase, and agreement with experiment confirms the pre- 
cision of the first principles methods. Finally, the transi- 
tion temperature of Polonium is calculated at pressure, 
comparing the Gibbs free energy of the different phases. 
Entropic and relativistic effects are delicately balanced 
and determine the stability of the a and f3 phases as a 
function of temperature. 

Density Functional Theory (DFT) [lj| and the per- 
turbation theory built on top of it (DFPT)[3] provide 
highly accurate calculations of the total energy and vibra- 
tional properties of materials from first principles. The 
implementation in the ABINIT package 15j is employed, 
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FIG. 1: The phonon band structure of SC polonium along 
high symmetry directions of the Brillouin Zone, including the 
spin-orbit interaction (solid line), and without it (dashed). 
Contrary to most normal metals, in Polonium spin-orbit de- 
creases nesting, removes Kohn anomalies, and hardens most 
phonon frequencies. 

with a plane wave basis (cut off at 435 eV) and pseudopo- 
tentials in the local density approximation [r|. Relativis- 
tic corrections are included through a pseudopotential of 
the HGH form[I3|. The corrections are decomposed into 
SO and SR contributions. The calculations confirm the 
interpretation of Ref. that the SR terms are sufficient 
to stabilize the simple cubic phase. Calculations without 
the SR relativistic terms will not be considered, because 
the of the structure of the pseudopotentials used, and 
because for heavy elements this is quite unrealistic. It is 
the spin-orbit term which has the largest influence on the 
shape of the Fermi surface (which is the central quantity 
affecting phonon structure) . The DFPT formalism with 
SO is described in Refs EU The thermodynamics 
of polonium are treated in the quasiharmonic approxi- 
mation, fitting volume and temperature dependencies as 
described by Grabowski et al.[20| The Brillouin Zone is 
sampled with 16 3 k-points for electrons and 8 3 q-points 
for phonons. 

The a phase of polonium is simple cubic, with a lat- 
tice constant of 3.345 A. The theoretical equilibrium lat- 
tice constant is 3.335 A, in good agreement, and will 
be used in the following. The phonon band structure 
of Po is shown in Fig. [TJ with and without the effect of 
the SO interaction. In a simple metal, the Fermi sur- 
face departs only slightly from a sphere. Kohn anoma- 
lies (KA)[21| appear for phonons with reciprocal space 
vectors q which are nested in the surface, i.e. those q 
which connect many pairs of points on the surface. These 
phonons are screened particularly strongly by the elec- 
trons, which renormalizes their frequency. Distorting a 
spherical FS, in particular because of the SO interaction, 
will always increase the nesting for certain q (sometimes 
dramatically as in Pb[l9|). In polonium, the opposite 
happens: due to the orbital nature and SC crystal struc- 



FIG. 2: Fermi surface of simple cubic Po, with (left) and 
without (right) the spin-orbit interaction. The SO interaction 
reduces the nesting of the Fermi surface, which in turn reduces 
the Kohn anomalies in Fig. [1] 



ture, the Fermi surface is more cubic than spherical (see 
Fig. [2]). As a result, the deformations due to the SO in- 
teraction will decrease the nesting, instead of increasing 
it. The small spherical hole pocket around T is also re- 
moved when the SO interaction is added. These types of 
effects had already been seen in other response proper- 
ties in Ref. [xoj | . The result is quite clear in the phonon 
bands of Fig. [T] when including the SO interaction (going 
from dashed to full lines), the dips due to Kohn anoma- 
lies are reduced, and the longitudinal (LA) frequency 
increases. The transverse (TA) frequencies increase or 
decrease in different regions of the Brillouin Zone, but 
also present KA which are reduced by the spin-orbit in- 
teraction. From the vibrational calculations for the a 
phase with and without SO, it is clear that inclusion of 
full relativistic effects is crucial to the correct represen- 
tation of the thermal properties of Po. On this basis, the 
thermodynamical behavior of a and /? Po (and their rela- 
tive stability) will be examined in the following including 
both SR and SO terms. 

Thermodynamical quantities are calculated from the 
phonon bands of the bulk metal. To capture effects of 
thermal expansion, full phonon calculations are also per- 
formed for a volume 3% larger than equilibrium, and 
the mode-dependent Griineisen parameters are extracted 
(see e.g. Ref. [2(] for standard formulae). The thermal 
expansion of Po is shown in Fig. [3] as a function of tem- 
perature; it is initially negative, which is found in many 
materials. For low temperatures the mode behavior near 
the Brillouin Zone center is crucial. The contribution 
from the smallest q vectors (up to 1/8 of Brillouin Zone) 
is therefore fit to explicit calculations at small q and in- 
tegrated numerically. 

The principal thermodynamical quantities for which 
there is experimental data are the thermal expansion co- 
efficient and the heat capacity. Brocklehurst et al.[22| 
found a thermal expansion coefficient of 23 • 10" 6 ± 
1.5 ■ 10~ 6 K _1 at 298 K. First principles calculations 
give 16 • 10~ 6 K _1 which is in reasonable agreement; the 
residual difference is probably due to anharmonic effects. 
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FIG. 3: Thermal expansion coefficient of Po as a function of 
temperature. Experiments^ give 23 ±1.5 10~ 6 K' 1 at 300 
Kelvin. 

The tabulated heat capacity of polonium at 300 K is 
Cp = 26.4 J mol -1 K -1 23], quite close to the ab initio 
value of 24.9 J mol -1 K _1 found here. The correction 
due to SO coupling is slight (+0.1 J mol -1 K _1 ). A 
fit to Murnaghan's [24j equation of state for a-Po gives 
a bulk modulus B = 39.45 GPa, volume 250 bohr 3 and 
B' = 4.89. 

We now turn to the thermodynamics of /3-Po, and the 
explanation of why this lower symmetry phase is favored 
at higher temperatures, examining the contributions to 
the free energy which stabilize a or (3 Po at different 
temperatures. The beta phase of Polonium is rhombo- 
hedral, with a lattice constant of 3.359 A and an an- 
gle of 98° 13'. It corresponds to a cube which has been 
slightly compressed along the (111) diagonal, and pos- 
sesses fairly low symmetry. As noted in previous stud- 
ies, the structure is intermediate between the SC and a 
body-centered-cubic phase, which is metastable at the 
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FIG. 4: The phonon band structure of (3 polonium along high 
symmetry directions of the Brillouin Zone, including the spin- 
orbit interaction 
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FIG. 5: Free energy differences (mHa) between j3 and a 
phases of Polonium, as a function of temperature (Kelvin). 
Positive values favor the stability of (3 Polonium. Top panel: 
total energy at K and thermally expanded volume. Middle 
panel: vibrational free energy at constant volume F v a. Bot- 
tom panel: volume-correction AF V ib (solid), electronic free 
energy at constant volume F e i (dash), and volume-correction 
AF e i (dot-dash). Differences are small below 600 K. It is the 
two vibrational free energy terms that determine the transi- 
tion temperature from a to f3 Po. 

equilibrium volume. For an angle of 60°, a face-centered- 
cubic phase is also metastable. For compressed volumes, 
calculations confirm the results of Ref. that a single 
rhombohcdral phase becomes stable for volumes smaller 
than about 0.93Vo (for a calculated pressure of 1 GPa), 
with an angle larger than 90°. Then, at compressions 
below 90% (for a calculated pressure between 1.5 and 
2 GPa), a second phase is also stabilized, with an an- 
gle smaller than 90° (as found in Ref. 0). Fitting an 
equation of state to the E(V) relation gives an analytical 
form for the basic variation of the energy with volume. 
This is not simple with /3-Po, for two reasons. First, 
phases with the equilibrium volume and unit cell angles 
between 90° and 100° are never even metastable with 
respect to a-Po. Second, fitting the equation of state 
to compressed rhombohedral phases is not useful, as an 
unrealistically large equilibrium volume is found. Given 
the good performance of the DFT methods for the SC 
phase, the experimental unit cell is used to compute the 
phonon bands of /3-Po, and the equation of state is fit for 
a constant rhombohedral angle (giving B = 40.15 GPa, 
volume 247 bohr 3 and B' = 4.49). Thermal expansion 
effects are then taken into account as for the a phase. 
This will be a very good approximation for volumes and 
temperatures close to our reference point, which is the 
experimental volume at around 348 K[6J. It has been 
checked that unstable modes do not appear in the vari- 
ous compressed or distorted structures. 

The phonon band structure of /3-Po is shown in Fig. 2J 
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along the same high-symmetry directions as in Fig. [T] (the 
SC nomenclature is used for simplicity) . The phonons are 
calculated at K (thermal expansion is later taken into 
account as for the a phase) , and would change somewhat 
due to anharmonic effects (see e.g. Ref . [iHh . Work in this 
direction is ongoing, but the current agreement with ex- 
periment suggests that anharmonic effects arc small. It 
has been checked that lowering the electronic tempera- 
ture does not lead to unstable phonons (in particular for 
the dip between L and T), such that the thermodynam- 
ics described by the K phonons is physical. The free 
energy is calculated for a and fl-Po, adding the ground 
state energy to the vibrational and electronic free en- 
ergies. The free energies are calculated at a second vol- 
ume V, and an additional correction is added, taking the 
thermal expansion into account (bottom panel of Fig. [5]): 
AF X (T) = {F x (V',T)-F x (Vo,T))-(V(T)-Vo)/(V'-Vo) 
where x is vibrational or electronic. The volume V{T) is 
obtained by inverting P — —dF/dV — 0, and inserted in 
the free energy to give the Gibbs free energy as a func- 
tion of temperature. The differences in the free energy 
components between the a and the (3 phases are shown 
in Fig. [5] In Ref. |9| the relativistic stabilization of ev-Po 
is calculated to be 0.15 mHa. Up to ^700 K, the to- 
tal energy difference between the phases (including the 
relativistic stabilization for both) is ca. -0.5 mHa. At 
350 K the vibrational entropy difference is of the order 
of 0.4 mHa, and the other terms also favor the (3 phase. 
The transition can thus be explained by the vibrational 
entropy terms in /3-Po overcoming the K total energy 
difference. 

Comparing the full Gibbs free energy of a-Po and (3- 
Po determines the transition temperature T c . The the- 
oretical T c (at pressure) from the a to the f3 phase is 
370 ±20 K. The error is estimated from the convergence 
of the free energy fit parameters, and does not take into 
account anharmonic corrections. This is in quite good 
agreement with experimental data which give a range 
of transition temperatures from 290 to 330 K. Omitting 
the thermal expansion or any of the free energy terms 
removes the transition and stabilizes a-Po for all tem- 
peratures. 

In summary, this Letter presents an ab initio study 
of the different thermodynamical phases of crystalline 
Po. The effect of spin-orbit coupling on the phonon band 
structure of simple cubic Po is first shown to be the op- 
posite of that on simple metals: as the initial Fermi sur- 
face (without SO) is strongly nested, the SO correction 
reduces the nesting and the resulting Kohn anomalies, 
and hardens most phonon frequencies. The inclusion of 
SO effects is essential to the correct representation of the 
phonons. Due to the delicate bonding characteristics of 
Po, whose low temperature simple cubic phase is stabi- 
lized by relativistic effects, increasing the temperature 
has the completely counter-intuitive effect of stabilizing 
a less symmetric phase (rhombohedral /3-Po). This tran- 



sition originates in the competition of relativistic stabiliz- 
ing effects for the a phase and slightly more favorable vi- 
brational entropy and thermal expansion for the (3 phase. 
A number of other elements are likely to have tempera- 
ture behaviors similar to Po: either due to relativistic 
or pressure effects, complex bonding environments open 
the way for radically different symmetry breaking as a 
function of temperature, when entropic effects are of the 
same order of magnitude as the energy differences be- 
tween bonding types. 
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